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ABSTRACT 

We analytically investigate the formation of an HII region in the accreting envelope 
of a newborn star. Special care is taken to examine the role of ionizing radiation force. 
This effect modifies velocity and density distributions and thereby affects the expansion 
of the HII region. As a result, the upper limit of the stellar mass imposed by the growth 
of an HII region around a forming star is increased by a larger factor than the previous 
estimate. In particular, for a star forming out of metal-free gas, this mechanism does 
not impose a firm upper limit on its mass. 

Subject headings: cosmology:theory — early universe — starsrformation — HII regions 


1. Introduction 

Stars are supposed to be born in the dense regions of gas clouds. The theory of the gravitational 
fragmentation of the gas clouds describes how these dense regions are formed inside the parental 
gas clouds. The typical mass scale of the dense region can be predicted by the Jeans mass of the 
cloud at the time of the fragmentation. However, a perturbation with a mass scale larger than this 
Jeans mass is always unstable. Therefore a dense region that is more massive than the typical Jeans 
mass can form depending on the initial fluctuation spectrum in the cloud. In this sense, the upper 
limit for the mass of stars cannot be obtained from the analysis of the gravitational fragmentation. 

Once, the origin of the observed upper limit of stellar mass (around 100M©) was attributed to 
the instability of massive stars due to the e-mechanism (e.g., Schwarzshild & Harm 1959). Since 
later studies revealed that this mechanism leads to moderate mass-loss rather than disruption of 
the stars (e.g., Appenzeller 1970), formation process, instead of the stability of massive stars, has 
been believed to limit the upper bound for the stellar mass (e.g., Nakano, Hasegawa, & Norman 
1995). 
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The star formation process is an accretion of ambient matter onto a protostar (stellar core) 
forming inside a protostellar cloud. In this scenario, the final mass of a star depends on how much 
mass the protostar can acquire. Although it is not yet clear what mechanism stops the accretion, 
increasingly strong feedback from a massive forming star is likely to halt the accretion. 

Larson & Starrfield (1971) were first to propose that a firm upper limit on the stellar mass is 
provided by the formation of an HII region around a forming star, as well as by the radiation force 
acting on dust grains in the accreting envelope. The formation of an HII region prevents further 
accretion in the following way: when the HII region reaches outer layers of a protostellar cloud, 
the temperature and then the pressure support surge by some large factor, so that the further 
infall of material is immediately halted. Larson & Starrfield (1971) concluded that both of those 
mechanisms set an upper limit of about 50 M 0 on the masses of Pop I stars forming out of the 
present interstellar medium. The radiation force onto dust grains was found by subsequent studies 
to set more stringent mass upper limit (Appenzeller & Tscharnuter 1974; Kahn 1974; Yorke & 
Kriigel 1977; Wolfire & Cassinelli 1987). Among them, Wolfire & Cassinelli (1987) suggested that 
even stars as massive as about 10 M 0 cannot form by accretion owing to this mechanism unless 
dust is significantly depleted. 

However, for so-called Pop III stars forming out of metal-free gas, this mechanism of the ra¬ 
diation force does not work because of the absence of dust grains. The e-mechanism drives even 
milder mass-loss for those stars than Pop I massive stars (Baraffe, Heger, & Woosley 2001). Ac¬ 
cording to recent simulations of the fragmentation of primordial clouds, the mass scale of fragments 
(protostellar clouds) is as large as 1O 3 M 0 (Bromrn, Coppi, & Larson 1999; Abel, Bryan, & Norman 
2000; Tsuribe 2001). Therefore, the growth of the HII region plays a crucial role in limiting the 
maximum mass of a Pop III star, if the upper limit by this mechanism is less than the mass scale 
of fragments. This motivates us to study here the upper mass limit of Pop III stars due to the 
formation of an HII region. 

In discussing the formation of an HII region, often the free-fall assumption has been imposed 
on the flow in the accreting envelope (e.g., Yorke 1986). The role of the momentum transfer to 
the gas due to ionizing radiation has been neglected. In the context of galaxy formation, however, 
its significance in dynamical and thermal evolution of the intergalactic medium has been pointed 
out by Haehnelt (1995). He showed, for the collapse of subgalactic clouds of < 1O 1O M 0 , even a 
radiation-pressure-driven bounce is possible. Taking this effect into account, in this paper, we study 
the formation of an HII region in the accreting envelope of a forming massive star. Although this 
effect is negligible as long as the radius of the HII region is small, it will turn out that inclusion 
of the ionizing radiation force alters the later evolution of the HII region. In the course of the 
expansion of the HII region, the flow becomes slower than the free-fall rate owing to the ionizing 
radiation force. The growth of the HII region is strongly suppressed by this effect. Consequently, 
the mass upper limit imposed by the expansion of the HII region is increased by a large factor. In 
particular, the formation of an HII region does not set any firm mass upper limit for Pop III stars. 
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In Section 2, we briefly summarize the relevant aspects of the formation of an HII region in a 
free-falling envelope around a newborn star. The effect of ionization radiation force is included in 
Section 3. Finally, we provide a summary and discussion in Section 4. 


2. Formation of an HII Region in a Free-Falling Envelope 


In this section, we review the expansion law of an HII region around a newly formed star 
(e.g., Yorke 1986) and discuss the resultant upper mass limit of a forming star assuming that the 
accreting flow is free-falling. 

Suppose an accreting flow onto a star of mass M * is taking place at an accretion rate M. The 
mass accretion rate M is determined by the temperature of the protostellar cloud and does not 
change very much except during the very late stage of the accretion phase where most of the mass 
of the envelope has accreted onto the central star. The relation between M and the sound speed 
c s in the protostellar cloud is given by 

M = C a cc |, (1) 

where C&cc is a non-dimensional number (see, e.g., Whitworth & Summers 1987). When the 
numerical value of M is needed, we use lO - 3 M 0 yr -1 corresponding to C acc = 1 and 400K which is 
a typical temperature of metal-free protostellar clouds (e.g., Yoneyama 1972; Carlberg 1981; Palla, 
Salpeter, & Stahler 1983). 

Here and throughout this paper, we assume that the accreting flow is spherically symmetric 
and steady. Then, the continuity equation yields 

M = —4:Trr 2 pu, (2) 


p and u are the density and velocity in the accreting flow at radius r. In addition, we assume in 
this section that the accreting envelope is free-falling: 


u = — 


2GM, 


Using equations (2) and (3), we obtain the density distribution 

M 


P = 


4ttV2GM* 


3/2 


( 3 ) 

( 4 ) 


In the HII region, we assume that the ionization and recombination are in equilibrium and that 
all the re-enritted ionization photons are consumed locally (a so-called on-the-spot approximation). 
Suppose Q ionizing photons are emitted per unit time from a sphere of radius i?; n surrounding 
the star. We regard the stellar surface at B, m as the base of the HII region for the moment (see 
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discussion in §4). For the sake of simplicity, we assume that the gas is purelycomposed of hydrogen 
and is fully ionized inside the HII region (i.e., the number density of H + n(H + ) = n(e) = p/m p ), 
hereafter. We obtain the radius Ru of the HII region from the ionization equilibrium: 


rR ii 

Q= an(R + )n(e)dV + 
J Rin 


M 
m v ’ 


( 5 ) 


where the recombination coefficient a = «b from the on-the-spot approximation. The first term 
on the right hand side means the recombination rate in the HII region. The second term means 
the neutral flux entering into the HII region. Note that no HII region forms if the second term 
on the right-hand side is larger than Q. In this paper, we neglect the M/m p term, assuming 
Q 3> M/m p = 3.8 x 10 46 m_3(s _1 ), where m_3 = M /10~ 3 Mqjt -1 . Using the relation (4), we 
obtain 

aM 2 

« = 5< 6 > 

Then, 

Ru = Rinexp(Q/Q^), (7) 

where 


Q 


(FF) 

crit 


aM 2 

87rG'm p 2 M* 

1 x 10 51 (Af*/100M Q )" 1 rh^ 3 (s^ 1 ). 


( 8 ) 

( 9 ) 


This shows that if Q exceeds Qc ri t , HII region expands exponentially. Soon after, the HII 
region reaches an outer layer whose gravitational binding energy is lower than the thermal energy 
of ionized matter. This results in the immediate expulsion of the outer layer and the halting of the 
accretion. 

/"pp'v 

When does the ionizing photon emissivity Q of the star reach the critical value Q crit • To be 
concrete, we consider here the case of Pop III stars. The ionizing photon emissivity 


Q ~ 1.6 x 10 50 (M*/1OOM 0 )(s _ 1 ) 


( 10 ) 


for very massive (M* > 300Mg) Pop III stars in the zero-age main sequence (ZAMS; Bromm, 
Kudritzki, & Loeb 2001). Comparing (10) with (9), we obtain the critical stellar mass above where 
Q exceeds 

M* ~3x lO 2 m_ 3 Af 0 . (11) 

The formation of a more massive star is increasingly difficult because of the growth of the HII 
region. 

Recall here that the exponential expansion law of the HII region (7), which plays a key role in 
determining the upper limit of stellar mass, is a direct consequence of the density distribution (4) 
and then the free-fall assumption (3). However, the gas is subject to momentum transfer due to 
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ionizing radiation and becomes slower than the free-fall rate within the HII region. This modifies 
the density law (4) and the expansion law of the HII region (7). Taking this effect into account, 
we study, in the next section, the growth of an HII region in the accreting envelope of a newborn 
massive star. 


3. Effect of Radiation Force on Expansion of an HII Region 


Hydrogen atoms in the HII region are ionized both by ionizing photons emitted directly by the 
star, and by re-emitted ionizing photons due to the recombination to the ground level. While the 
former is radially outward-directed, the re-emitted photons become almost isotropic and contribute 
little to the radiation force in the radial direction. For the calculation of radiation force, we must 
subtract the contribution from the ionization by locally emitted photons from the total ionization 
rate f3n( H). The rate of the former equals the recombination rate to the ground level aqn(H + )n(e) 
(per unit volume) according to the on-the-spot approximation, while the latter equals the total 
recombination rate ap L n{B + )n{e) by ionization equilibrium. Hence, the radiation force by ionization 
radiation per unit mass is given by 


/ion = —^-an(H + )n(e) / p, 
c 


( 12 ) 


where, a = «b again. Since we are 
energy of ionizing photons is hv i on ~ 
radius 


first introduced by Haehnelt (1995), 
ionizing radiation force and gravity. 4 


considering stellar UV sources of thermal spectra, the mean 
hv o = 13.6eV. For later convenience, we define the Haehnelt’s 


i?H = 


ahv\, 


4:TTcGrn'p 


26pc, 


(13) 


which is a fundamental length scale for the balance between 
Finally, the expression of the ionizing radiation force becomes 


/ion = 4 ttGRhP. (14) 

Note that we have neglected other photon extinction processes such as electron scattering, pho¬ 
todissociation of H“, etc. This assumption is valid as long as the HII region outside R\ n is optically 
thin to those extinction processes. 

We consider here an HII region with the steady radius Ru determined by the ionization equilib¬ 
rium (5) as in the previous section. The I-front moves outward slowly as the stellar mass, and then 
the ionization photon emissivity from the star, increase. If the flow velocity across the ionization 
front (I-front) is more than twice the sound speed of ionized matter, the I-front becomes a (weak) 


4 Comparing radiation force /i 0n = ‘inGRHp from a point source with gravity / grav = G(4npr 3 /3)/r 2 exerted 
by homogeneously distributed matter with density p around the source, Haehnelt (1995) found the radiation force 
dominates the gravity at r < 3 Rh- 
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R-type (Kahn 1954). Since we consider the steady radius of the HII region, the velocity across the 
I-front equals that of the flow at the edge of the HII region. As long as the HII region is not large, 
the flow velocity is supersonic and fast enough to keep the I-front an R-type. In fact, as we will 
show later, the halting of accretion occurs while the I-front is still R-type. Hence, we treat here 
only the HII region bounded by R-front and do not consider later evolution. 

Including the ionization radiation force, the momentum equation for the steady flow yields 

du GAR 

u— =-tt -1- 47rGi?H/9 (15) 

dr r z 

in the HII region ( R\ n < r < An)- The thermal pressure term can be neglected since we consider 
only supersonic flow. 


Using the density-velocity relation in the steady accretion 


M ! 
P = ^- U) ’ 


(16) 


and introducing the kinetic energy density of the flow K = u 2 / 2, we obtain from equation (15) 


dK = _GAR ( _ [K~^\ 

dr r 2 ( V K J 


(17) 


Here we define the critical kinetic energy A cr it and the critical flow velocity n cr it ; 

1 1 

Acrit = 2 ^crit = ^(-^H/^acc) > 


and the accretion timescale 

face = AA/M. 


(18) 

(19) 


Given the value of K(= An) at the outer boundary An of the HII region, we can solve the 
momentum equation (17) inward. Depending on whether the initial value A"n is larger than, equal 
to, or smaller than the critical value A' cr it; solutions of the equation (17) are classified into three 
types. In Figure 1, the solutions for An/A cr ; t = 1.1,1.01,1,0.99, and 0.9 are shown. If An is 
larger than the critical value K C rit, the right hand side of (17) is always negative. Then A becomes 
larger and larger as the gas flows inward: The flow is accelerated. We call these types of solutions 
accelerating solutions. If An = A cr it, A remains always at this value (the uniform velocity solution), 

u(r) = Merit = (20) 

^acc 


In other words, the ionizing radiation force exactly balances the gravity everywhere in the HII 
region. From equations (2), (18) and (19), the density distribution of the uniform velocity solution 
is given by 


P = 


M* 

47tAh 


r 


-2 


( 21 ) 
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Similarly, if An < A cr i t , the flow continues to decelerate (decelerating solutions). As can be seen 
in Figure 1, the deceleration is so rapid that the accretion is halted by the radiation force without 
reaching the star shortly after the critical flow is attained. 


As the outer boundary condition, we assume that the flow is free-fall outside the ionizing 
region. We further assume that the flow velocity is continuous at the I-front, since velocity jump 
at the R-front is small. Then, 


An = 


GM* 

An 


( 22 ) 


the kinetic energy An at the outer boundary of the HII region is related to the radius An of the 
HII region. 


The flow equation (17) can be solved analytically. The solution satisfying the boundary con¬ 
dition (22) is 

— = A + 2^K~ t (VK - /An) + 2A crit ln( (if A n + A crit ), (23) 

r v/An - VAcnt 

and 

A = A crit (if A„ = A crit ). (24) 

In the limit of small radiation force (i.e., A cr ;t —► 0), only the first term on the right hand side of 
equation (23) remains and the flow becomes free-fall. The second and third terms are corrections 
due to the radiation force. 


From the above characteristics of the flow equation (17), the growth of the HII region can be 
described as follows: As long as the the forming star is not massive enough to emit a large quantity 
of ionization photons, the HII region around it is compact. Since the radius An is small, An > A cr ;t 
in this stage. The flow inside the HII region is hardly affected by the ionizing radiation force. As 
the star grows in mass and begins to emit a large quantity of ionizing photons, the radius of the 
HII region becomes larger and larger. When the radius of the HII region reaches the critical value 

2 GM 

Ancrit = " 77 / = 0.59 x lO 3 A 0 (M*/lOOM Q ) 3 m/, (25) 

V-^H/^acc ) 

a uniform velocity flow is attained inside the HII region. Further growth of the HII region would 
result in non-steady accretion and the accretion would eventually halt. We cannot discuss the later 
evolution of the accreting flow within the framework of steady accretion. According to hydrody- 
namical calculations for the accretion being halted by the radiation force onto dust, the accretion 
rate begins to decrease with some oscillations when the luminosity reaches the critical value (Yorke 
& Kriigel 1977). We speculate that a similar phenomenon would occur also in our case. 

We justify here the assumption that the I-front remains in R-type until it reaches the critical 
radius An cr it. When the radius of the HII region reaches An cr it, the flow velocity in the HII region 
is Ucrit (see eq. 20). The I-front is R-type if 


l^critl > 2c s n 


( 26 ) 



where c s n is the sound speed of the ionized matter. Then, our assumption of R-front is justified as 
long as 

tacc <^^2X 10 6 yr, (27) 

2c s n 

where for Cgii we used the value at 10 4 K. The right hand side of equation (27) is as large as the 
lifetime of a massive star. Recall that the accretion time t acc means the time needed to form a 
star of mass M* at an accretion rate M. Since the halting must occur within the stellar lifetime, 
the condition (27) is always met. Therefore, in our case, we can conclude that the I-front remains 
R-type until the halting of the accretion. If the mass accretion rate decreases significantly, the 
formally defined accretion timescale f acc (eq. [19]) can be larger than the stellar lifetime, and hence, 
the condition (27) would be violated within the stellar lifetime. However this decrease of the mass 
accretion rate is expected only in the very late stage of the mass accretion phase where the mass 
of the envelope is significantly depleted. Therefore the mass of the star eventually becomes on the 
order of the mass of the parental cloud. That is, the halting mechanism due to the ionizing photons 
considered here is not important in this case. 


Next, we estimate the ionizing photon emissivity needed for attainment of the uniform velocity 
solution. Substituting the density distribution at the critical flow (21) into the equation of ionization 
equilibrium (5), we obtain the critical ionization photon emissivity 


Q crit — 


olM 2 


1 


87rm|A" crit R in 
fhv Km GMl 


Rllcrit 

1 1 


i?H -Rin Rllcrit 


(28) 

(29) 


above in which the flow is decelerated in the HII region. We have neglected the accretion term 
M /mp as in §2. 


In equation (28), if the right hand side is dominated by the first term (i.e., R in <C Riicrit), 


Qcrit = (—r^ = 0.64 x lO 53 (R in /lOR 0 )- 1 (M,/lOOM 0 ) 2 . (30) 

C -*Mn 

As a typical value of i?i n , we take 10 -Rq, which is approximately the radius of ZAMS Pop III stars 
of 5OOAf 0 (El Eid, Fricke, & Ober 1983). 

Although we describe both Q cr it (eq. [28]) and Q) rit ; in §2 (eq. [9]) as the critical ionizing 
photon emissivities, mechanisms of halting the accretion are different. In §2, the accretion is 
halted by the gas pressure in the HII region, whereas here what halts the accretion is the radiation 
pressure of ionizing photons. The inclusion of ionization radiation force results in enhancement of 
<2 CT it because of an increase in the density (and hence, the recombination rate) due to the reduced 
flow velocity. 


So far, we have only considered ionizing radiation as a source of radiation force. However, for 
massive stars the luminosity approaches the Eddington limit, and then the radiation force due to 
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electron scattering becomes important. We can include this effect in our scheme without difficulty. 
If we use G e s = G( 1 — T) and Rn,eS = Rh/0- ~ T) (where T = L/L^d) instead of G and i?Hi 
the flow equation including the radiation force due to electron scattering retains the same form 
as equations (15) and (17). Note that G in (22) is not replaced since the gravity outside the HII 
region is not altered. After the same procedure as that leading to (28), the critical emissivity of 
ionization photons becomes 


Qcrit 


hu ion i GM* (l-r ) 2 

C Ru R'm 


Z?IIcrit 


(31) 


The critical ionizing photon emissivities Q C rit given by equations (28) and (31), namely the 
critical values with and without the effect of the electron scattering, are both depicted in Figure 
2 as functions of the central stellar mass. In this figure, we have used M = lO^Mgyr^ 1 and 
R[ n = 10-Rq as before. Also shown are the ionizing photon emissivity Q of ZAMS Pop III stars 
and the critical value Q { cvit obtained by assuming the free-fall (eq. [9]). While Q exceeds Q(. rit 
at about 300 Mq as already mentioned in §2, Q cr it is about two orders of magnitude larger than 
the actual value Q even at 10 3 M Q for our assumed value of R,\ n = 10 Rq, even though the electron 
scattering reduces Qcrit considerably. Recalling that Q C rit oc , we can conclude that Q C rit > Q 
for M* < 1O 3 M 0 unless Ri n becomes as large as lO 3 iZ 0 and Q remains the same value. We consider 
that this situation does not occur in reality as long as we regard R.\ n as the stellar surface. Therefore, 
the growth of a HII region does not prevent the formation of such massive stars as 1O 3 M 0 , 


4. Summary and Discussion 

To estimate the mass upper limit of metal-free stars imposed by the formation of an HII region, 
we have found solutions for spherically symmetric and steady accretion flow onto a star emitting 
ionizing radiation. 

The behavior of the flow is determined by the velocity of the flow entering the HII region. If the 
flow velocity at the edge of the HII region goes beyond (below) a critical value, the flow is always 
accelerating (decelerating, respectively) in the HII region. When this value is equal to the critical 
value, the velocity remains constant in the HII region. In this critical flow, the radiation force 
due to ionizing photons exactly balances the gravity. We applied those solutions for a compact 
HII region forming around an accreting star. As the HII region grows in radius, the accreting 
flow evolves from an accelerating solution to a critical one. Soon after the critical flow is reached, 
accretion is halted by radiation force due to the ionizing radiation. However, even stars as massive 
as 10 3 AZq are unable to emit ionizing photons sufficient to halt the accretion. The halting by gas 
pressure is even more difficult. Therefore, contrary to the previous expectation, the formation of 
the HII region does not impose a stringent mass upper limit (at least up to 1O 3 M 0 ) on metal-free 
stars. 

We have estimated the upper limit of the mass of Pop III stars in relation to the formation of 
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an HII region. Note that other effects (e.g., stellar wind, mass outflow due to pulsational instability, 
etc.) could grow in importance and might decrease the upper mass limit of Pop III stars below the 
value obtained in this paper. We leave detailed studies on these topics for future work. 

Also, in this paper, we have adopted very simplified assumptions, for example spherical sym¬ 
metry, steady accretion, etc. Here, we discuss other complexities and possible deviations from our 
picture. 

We have considered only photoionization and electron scattering as sources of radiation force. 
Here, we mention briefly the radiation force due to Lya photons. Lya photons are emitted from 
recombination in the HII region. The flow inside the HII region is exerted by the radiation force 
due to Lya photons. We can extend our theory easily to include the Lya pressure , since, in the HII 
region, the radiation force due to the Lya photons is proportional to that due to ionization photons 
(equation 12 of Haehnelt 1995; see also Braun & Dekel 1989; Bithell 1990). However, according to 
Haehnelt(1995), the radiation force due to Lya photons is at most of marginal importance relative 
to that due to the ionizing photons. Thus, for the sake of simplicity, we have not included it. The 
HII region is surrounded by an HI layer, which is very optically thick to Lya photons. Without 
dust grains, which absorb Lya photons and reemit them as infrared photons, the Lya photons 
must diffuse out through the HI layer. In the course of this, the HI layer is pushed outward by 
those photons. Doroshkevich & Kolesnik (1976) argued that this mechanism expels the HI layer 
soon after the HII region is formed and thereby limits the mass of stars below 10 Mq. However, 
Harrington (1973) showed that, even without dust, the two-photon emission process decreases the 
Lya photon density drastically, and the Lya radiation force is not dynamically important in an 
HI layer surrounding an HII region. In our case, supersonic motion in the accreting envelope also 
decreases the Lya photon density. Considering these facts, it is likely that the Lya radiation force 
does not play a significant role in our case. Therefore, we chose to neglect it here and assumed the 
free-fall outside the HII region. 

We have identified the base of an HII region as the stellar surface, and have taken the typical 
value of the stellar radius lOi?© as the inner boundary radius 7?; n of the HII region. However, when 
the accretion rate is high, a photosphere may be formed in the accreting flow (Wolfire & Cassinelli 
1986; Stahler, Palla, & Salpeter 1986). In this case, R m should be taken as the radius of the 
photosphere, since we assume the HII region is optically thin to continuum absorption except for 
photoionization outside R\ n . In spite of higher R \ n , this effect generally works towards the smaller 
HII region, since Q drops (Wolfire & Cassinelli 1986). Hence, our conclusion of higher upper mass 
limit than the former estimate remains the same. 

The spherical symmetry is clearly an oversimplification if the accreting matter has large angular 
momentum and an accretion disk is formed. Even in this case, the spherical symmetric flow is a 
good approximation outside the centrifugal radius. Suppose here that inside the centrifugal radius, 
there is a “cavity” around the disk. In this case, /?j n should be taken as the radius of the cavity. If 
Q is unattenuated inside the cavity, the halting of the flow outside becomes easier because of higher 
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i?i n and then lower Q C rit (see eq. [28]). The accretion through the disk might continue, however 
(Nakano 1989). Those issues are still too speculative and beyond the scope of this paper. 

Although the accelerating solution is Rayleigh-Taylor stable, it might be unstable to pertur¬ 
bations: if a portion of the flow becomes slightly slower than the rest, it becomes slower and slower 
relative to the average flow owing to the increased radiation force. In this case, density irregularities 
or blobs could be formed in the flow, and our spherically symmetric solution might be regarded 
as an approximate description of the average flow. Further study of the stability and dynamical 
evolution of the flow will be interesting. 


We thank H.Sato for helpful comments and V.Bromm for providing the data of Pop III ZAMS. 
This work is supported in part by Research Fellowships of the Japan Society for the Promotion of 
Science for Young Scientists, grant 6819. 


REFERENCES 

Abel, T., Bryan, G. L., & Norman, M. L. 2000, ApJ, 540, 39 

Appenzeller, I. 1970, A&A, 9, 216 

Appenzeller, I., & Tscharnuter, W. 1974, A&A, 30, 423 

Baraffe, I., Heger, A., & Woosley, S. E. 2001, ApJ, 550, 890 

Bithell, M. 1990, MNRAS, 244, 738 

Braun, E., & Dekel, A. 1989, ApJ, 345, 31 

Bromrn, V., Coppi, P. S., & Larson, R. B. 1999, ApJ, 527, L5 

Bromrn, V., Kudritzki, R. P., & Loeb, A. 2001, ApJ, 552, 464 

Carlberg, R. G. 1981, MNRAS, 197, 1021 

Doroshkevich, A. G., & Kolesnik, I. G. 1976, Soviet Astron., 20, 4 
El Eid, M. F., Fricke, K. J., & Ober, W. 1983, A&A, 119, 54 
Haehnelt, M. G. 1995, MNRAS, 273, 249 

Harrington, J. P. 1973, Memoires Societe Royale des Sciences de Liege, 6, 437 
Kahn, F. D. 1974, A&A, 37, 149 

Kahn, F. D. 1954, Bull. Astron. Inst. Netherlands, 12, 187 



12 - 


Larson, R. B., & Starrfield, S. 1971, A&A, 13, 190 
Nakano, T. 1989, ApJ, 345, 464 

Nakano, T., Hasegawa, T., & Norman, C. 1995, ApJ, 450, 183 

Osterbrock, D. E. 1986, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei (Mill Valley: 
University Science Books) 

Palla, F., Salpeter, E. E., & Stahler, S. W. 1983, ApJ, 271, 632 
Schwarzschild, M., & Harm, R. 1959, ApJ, 129, 637 
Stahler, S. W., Palla, F., & Salpeter, E. E. 1986, ApJ, 302, 590 

Tsuribe, T. 2001 in ASP Conf. Ser. 222, The Physics of Galaxy Formation, ed. M. Umemura, & 
H. Susa (San Francisco:ASP), 27 

Whitworth, A. & Summers, D. 1985, MNRAS, 214, 1 

Wolfire, M. G., & Cassinelli, J. P. 1986, ApJ, 310, 207 

Wolfire, M. G., & Cassinelli, J. P. 1987, ApJ, 319, 850 

Yoneyama, T. 1972, PASJ, 24, 87 

Yorke, H. W. 1986, ARA&A, 24, 49 

Yorke, H. W., & Kriigel, E. 1977, A&A, 54, 183 


This preprint was prepared with the AAS IAILX macros v5.0. 




- 13 - 




r /-^II,crit r /^II,crit 


Fig. 1.— Solutions of accretion flows in an HII region, (a) the kinetic energy K and (b) density p 
distributions for flows with Ku/K cr i t =1.1 (dash-dotted), l.Ol(long-dashed), 1 (solid), 0.99(dashed), 
and 0.9 (dotted). The thick solid line depicts the free-falling flow. The density normalizer pncrit is 
the density of the uniform-velocity solution (i.e., Ku/K cr it = 1) at r = R\\. The other normalizers 
can be found in the text. 










- 14 - 



Fig. 2.— The ionizating photon emissivity of Pop III stars Q vs the critical value for halting of 
the accretion Q C rit as a function of the stellar mass. The solid line shows the actual ionization 
photon emissivity of ZAMS Pop III stars. The critical values Q c .rit are shown for both cases with 
(long-dashed; eq. [31]) and without (short-dashed; eq. [28]) the effect of electron scattering. Also 
shown is the critical value Q rrn under the free-falling assumption (dotted; eq. [9]). We assume 
the following values of parameters: M = 10 -3 MQyr -3 and R m = 10-Rq. The values of Q and T 
are adopted from Bithell (1990) for M* < 500M* and Bromm et al.(2001) for higher masses. 





